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Decays of the T(5S) resonance to channels with B + and B° mesons are studied using a 23.6 fb -1 
data sample collected with the Belle detector at the KEKB asymmetric-energy e + e~ collider. Fully 
reconstructed B + J/ipK + , B° -> J/tpK* , B + -> D°ty + and B° D~tt + decays are used to 
obtain the charged and neutral B production rates per bb event, f(B + ) = (72.llg'g ± 5.0)% and 
f(B°) = (77.0lg;6 ± 6.1)%. Assuming equal rates to B + and B° mesons in all channels produced 
at the T(5S) energy, we measure the fractions for transitions to two-body and three-body channels 
with B meson pairs, /(BB) = (5.5 ± 0.4) %, f(BB* + B* B) = (13.7 ± 1.3 ± 1.1) %, f{B*B*) = 
(37.5 ±\\ ± 3.0) %, f{BB 7r) = (0.0 ± 1.2 ± 0.3)%, f(BB*ir + B* Bit) = (7.3 1\\ ± 0.8) %, and 
f(B*B*ir) = (1.0li;3±0.4)%. The latter three fractions are obtained assuming isospin conservation. 

PACS numbers: 13.25.Gv, 13.25. Hw, 14.40.Pq, 14.40.Nd 



New aspects of beauty dynamics can be explored us- 
ing the large data sample recently collected by Belle at 
the center-of-mass (CM) energy of the T(5S) resonance 
(also referred to as T(10860)). Decays of bottomonium 
states with masses higher than the mass of the T(4S) 
are almost unexplored with only a very limited num- 
ber of experimental and theoretical studies to date. At 
the T(5S) energy a bb quark pair can be produced and 
hadronize into various final states, which can be classi- 
fied as two-body B° s (B° s B° L B s Bl, B*B°, B*B*) [IS], 
two-body B(BB, BB*, B*B, B*B*), three-body (BBtt, 
BB*n, B*Bn, B*B*n), and four-body (BBitit) chan- 
nels. Here B denotes a B + or B° meson and B denotes 
a B~ or B° meson. The excited states decay to their 
ground states via B* — > Bj and B* — > B®j. Moreover, a 
bb quark pair can also hadronize to a bottomonium state 
accompanied by 7r, K or r\ mesons, for example through a 
T(5S)— > Y(lS)7r + 7r~ decay In addition, initial-state 
radiation (ISR) can affect the final states listed above 
and must be taken into account Fractions for all 
of these channels provide important information about 
fe-quark dynamics. 

The first study of B ^production at the T(5S) was per- 
formed by CLEO [H, |(| using a 0.42 fb -1 data sample. 
They found the fraction of events with B + /° pairs to be 
(58.9 ± 10.0 ± 9.2)%. CLEO interpreted the remaining 
(41.1 ± 10.0 ± 9.2)% as the fraction of events with B" 
mesons. Within the large uncertainties this fraction is 
1.6cr larger than the B° s event fraction f s = (18.0 ± 1.3 ± 
3.2)%, directly measured by Belle [2j . Both values for this 
fraction were obtained assuming that contributions from 
channels with bottomonium states are negligibly small. 
Among the _B+/° pair events CLEO found the two-body 



fractions for the B*B* and BB* + B*B channels to be 
(74 ± 15 ± 8)% and (24 ± 9 ± 3)% @, respectively. The 
BB channel and multibody channels were not observed 
and corresponding upper limits were set. 

Several theoretical papers have been devoted to T(5S) 
decays to final states with two-body B® and B + /° pairs 
[7HlO|. The B*B* channel is predicted to be dominant 
with the fraction over all bb events within the range 
(30-69)% [l[l0j]. In these model calculations the other 
two possible channels have smaller fractions with predic- 
tions covering a broad range. Multibody channels have 
also been theoretically studied [ll|, [l2[; the three-body 
fractions are found to be about two or three orders of 
magnitude smaller than the two-body fractions. Inter- 
esting information about a possible gluonic component 
of the T(5S) can also be obtained from measurements of 
the three-body decays [l3| , if 200 or more events can be 
reconstructed in a three-body channel. 

Here we fully reconstruct the decay modes B + — > 
J/ipK+, B a J/i>K*°, B+ -> D a ir+ (using two D° 
modes), and B° — > D~ir + . Charge-conjugate modes 
are implicitly included throughout this work. The re- 
constructed modes have large and precisely measured 
branching fractions LJ| and contain only charged par- 
ticles. 

The data were collected with the Belle detector [l5| at 
KEKB (l6j . an asymmetric-energy double storage ring 
e + e~ collider. This analysis is based on a sample of 
23.6 fb" 1 taken at the T(5S) CM energy of -10867 MeV 
and containing A^ (5S) = (7.13 ± 0.34) x 10 6 produced 
events @- The Belle detector is a general-purpose mag- 
netic spectrometer described in detail elsewhere [lj| . 

Charged tracks are assigned as pions or kaons based on 
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a likelihood ratio C^/^ = Ck/ {£k +C w ), which includes 
information obtained from the three Belle particle identi- 
fication detector subsystems [3]. The identification effi- 
ciency for particles used in this analysis varies from 85% 
to 92% (91% to 98%) for kaons (pions). The electron and 
muon identification requirements are described in Refs. 
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The K*°, D°, D 
structed in the K*° 

K + TT 

J/t/j 



fl+fj, 



and J/ip candidates are recon- 

-> K+TT", D° -> K+TT', D° -> 

D~ — > A" + 7r~7r~, J/ip — > e + e~ and 
modes. We require the invariant masses 



to be in the following intervals around the nominal 
masses: ±50MeV/c 2 for K*°, ±10MeV/c 2 for D° and 



D~, ±30MeV/c 2 for J/ip -)• 



andi?° MeV/c 2 for 



J/ip — > e + e~. A vertex- and mass-constrained fit is ap- 
plied to J/tp, D° and D~ candidates to improve the B 
signal resolution. 

B decays are fully reconstructed and identified us- 
ing two variables: the energy difference AE = E™ — 
£™ m and the beam-energy-constrained mass Afbc = 

^(^bcam) 2 - (Pb M ) 2 ' whcrc E b M and P™ arc the en- 
ergy and momentum of the B candidate in the e + e~ 
CM system, and -E^L is the CM beam energy. The in- 
termediate two-body and multibody channels with £> + /° 
pairs cluster in distinct regions of the Mh c and AE 
plane. However, all channels are distributed along a 
straight line described approximately by the function 
AE = tub — Mbo where raj is the nominal B mass 
(fixed to 5.28 GeV/c 2 ). 

After all selections the dominant background is from 
e + e~ — > qq continuum events (q = u,d,s, or c). Events 
with B mesons tend to be spherical, whereas continuum 
events are expected to be jet-like. To suppress contin- 
uum background, we apply topological cuts. The ratio 
of the second to the zeroth Fox- Wolfram moments [l9| is 
required to be less than 0.5 for the low background final 
states with a J/ip, and less than 0.4 for all others. The 
angle in the CM between the thrust axis of the particles 
forming the B candidate and the thrust axis of all other 
particles in the event, ^ hr , must satisfy |cos(9* hr | < 0.9 
for the final states with a J/ip, and |cos0j hr | < 0.75 for 
all others. More than one B + /° candidate per event is 
allowed, however the probability of multiple candidates 
is less than 1% for the modes used here, where the final 
states contain charged particles only. 

The two-dimensional M^ c and AE scatter plots for 
the B+ J/ipK+, B° J/V>A*°, B+ -4 D°ir+ 



D- 



(L>° -> K+tt- and L>° — » K+tt+tt-tt-) and B 
modes are obtained and are shown in Fig. 1. Events 
are clearly concentrated along the line AE = mj — A/b c 
corresponding to B pair production. The signal regions 
shown in Fig. 1 are restricted to a ±30 MeV interval in 
AE (corresponding to (2.5-4.0) a for the studied modes) 
and to the kinematically allowed 5.268 GcV/c 2 < M bc < 
5.440 GeV/c 2 range. The location of specific channels 



within the signal bands will be discussed below and is 
shown in Fig. 2a. 
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FIG. 1: The M bc and AE scatter plots for (a) B + -»■ 
J/ipK+, (b)_B° -> J/1PK* , (c) B+ -> D°(K + TT-)ir+, 
(d) B + -> D°{K + tv + it-iv-)it + , and (e) B° ->■ D--K+ de- 
cays. The bands indicate the signal regions corresponding 
to the intervals 5.268 GeV/c 2 < M bc < 5.440 GeV/c 2 and 
\AE + M bc -m B \ < 0.03 GeV. 

We use inclined AE + Mb c — ms projections of the 
two-dimensional scatter plots for all events within the 
range 5.268 GcV/c 2 < M hc < 5.440 GeV/c 2 to obtain 
integrated B decay event yields. We fit these distribu- 
tions with a function including two terms: a Gaussian 
to describe the signal and a first-order polynomial to de- 
scribe background. The regions where cross-channel B 
backgrounds can contribute are not used in the fits. 

Using the fit results, the charged and neutral B pro- 
duction rates per bb event are obtained from the formula: 



f(B^) = Y^ x I (JVj< 5S > x e B -x x B B -+x) , (1) 

where Yg\ x is the event yield obtained from the fit for 
a specific mode B — > X, es^x is the reconstruction ef- 
ficiency including intermediate branching fractions, and 
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Bb^x is the corresponding B decay branching fraction 
14 1. The event yields, efficiencies, and production rates 
are listed in Table I. 

TABLE I: Event yields obtained from fits, efficiencies, and B 
production rates f(B + ^°). Efficiencies include intermediate 
J ftp, K*° , D° and D~ branching fractions. 



Decay mode 


Yield 


Efficiency, % 


f{B+>°),% 


B + -> J/ipK + 


221±il 


3.41 


89.0t6: : i ± 8.0 


B° ->■ J/^jK* 


105 ± ff 


1.30 


85.3±|;| ±8.8 


B+ D°(Ktt)it + 


215 ±21 


0.97 


64.0 ±6.2 ±4.9 


B+ -> D°(K3tt)tv + 


275 ± 32 


1.17 


68.3ig;" ±6.4 


B° -»• + 


247 ± 25 


1.80 


72.9 ±7.4 ±6.4 



The systematic uncertainties include those due to the 
determination of the number of bb events (4.7%), charged 
track reconstruction efficiency (1% per track), particle 
identification (0.5-1.0% per it and K, 2% per electron, 
3% per muon), J/ip, D° and D~ mass cut efficiencies 
(2%), signal and background modeling in the fit proce- 
dure (2%), MC statistics in efficiency determination (1- 
2%), shape of the B meson angular distribution relative 
to the beam axis direction in the CM system (1%), and 
PDG branching fractions (3-5%). All systematic uncer- 
tainties are combined in quadrature to obtain the total 
systematic uncertainty. The same set of uncertainties 
is used below to obtain the total systematic uncertain- 
ties for values averaged over several B modes, however 
the correlated and uncorrelated uncertainties are treated 
separately in this case. All uncorrelated uncertainties 
from all B modes are varied individually to obtain their 
contribution to the total systematic uncertainty of the 
averaged value, while correlated uncertainties are varied 
simultaneously for all channels. 

Using the rates shown in Table 1, the average pro- 
duction rates f(B+) = (72.1±|;| ± 5.0)% and f(B°) = 
(77.0ig;g ± 6.1)% are obtained. As expected, these rates 
are equal within uncertainties. The average of charged 
and neutral B modes is (73.7 ± 3.2 ± 5.1)%. Within un- 
certainties this rate is consistent with the CLEO value of 
(58.9 ± 10.0 ±9.2)% @. 

Since the f(B + ) and f(B°) values are consistent with 
isospin symmetry, it is reasonable to assume that equal 
numbers of B + and B mesons are produced in all pos- 
sible channels. Therefore the five B decay modes are 
treated simultaneously everywhere below. First, events 
in the signal bands shown in Fig. 1 are projected onto 
Afbc- To describe combinatorial backgrounds under the 
signal, we define sidebands, which have the same shape as 
the signal region, but are shifted by ±70 MeV in AE, and 
project similarly onto M\, c . Assuming that background is 
distributed linearly in AE, we model the combinatorial 
background by taking the average of the two sideband 
Mbc distributions. 

Figure 2(a) shows the Mb c distributions for the pos- 



sible two-, three- and four-body channels obtained from 
MC simulation for the B° — > D~ir + decay. The three- 
body and four-body channels are simulated assuming 
a pure phase-space decay model. Figure 2(a) shows 
that the two-body channels are well separated. The 
three-body BBn, BB* tt + B*Btt and B*B*ir chan- 
nels have broader distributions in the higher mass region 
M bc >5.35 GeV/c 2 . The four-body BBtttt channel is 
the cross-hatched peak on the right side of Fig. 2(a). 




M bc (GeV/c') 



M bc (GeV/c' 



FIG. 2: (a) MC simulated M hc _ distributions for the B° -> 
D~tt + decay for BB, BB* + B*B, B*B* and BBtttt channels 
(cross-hatched histograms from left to right), and also for the 
three-body channels BB* tt + B*Btv (plain histogram), BBir 
(dotted) and B* B* tt (dashed). The distributions are normal- 
ized to unity, (b) M^ c distribution in data after background 
subtraction. The sum of the five studied B decays (points 
with error bars) and results of the fit (histogram) used to ex- 
tract the two-body channel fractions are shown. The fitting 
procedure is described in the text. 

The matrix elements responsible for the three- and 
four-body decays are not known, and the rates of the 
three- and four-body contributions cannot be obtained 
in a model-independent way from a fit to these Mb c 
distributions. We therefore restrict the fit to the re- 
gion 5.268 GcV/c 2 < M bc < 5.348 GcV/c 2 to ex- 
tract the two-body channel fractions. To obtain the 
sum of the other contributions in the large Mb c re- 
gion, only the events from the interval 5.348 GcV/c 2 < 
Mbc < 5.440 GeV/c 2 arc selected. For these distribu- 
tions we apply a fit procedure similar to that used above 
to obtain the full f(B + /°) production rate for the entire 
interval 5.268 GeV/c 2 < M bc < 5.440 GeV/c 2 . The to- 
tal signal yield, summed over all modes and channels in 
this region, is 228.7 j^J. 

The fitting procedure used to extract the production 
rates of the two-body decay channels treats the five signal 
and five sideband Mb c distributions simultaneously. The 
following four components with fixed shapes and floating 
normalizations are included in the fit for each distribu- 
tion: BB, BB* + B*B, B*B* , and combinatorial back- 
ground. The shapes of the signal components are taken 
from MC simulation, and those of combinatorial back- 
ground are modeled by the sideband data. The normal- 
ization parameters (channel fractions) for the two-body 
components arc constrained to be equal for all five stud- 
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ied B decays. The sum of background-subtracted Mb c 
distributions for the five studied B decays is shown in 
Fig. 2(b), where the result of the likelihood fit in the 
region 5.268 GeV/c 2 < M hc <5.348 GcV/c 2 used to ob- 
tain two-body channel fractions is superimposed. The 
fractions obtained from the fit are listed in Table II. The 
systematic uncertainties include all those described above 
as well as uncertainties due to the Mh c signal shape mod- 
cling (3%). 

TABLE II: The fractions of two-body channels with B + and 
B mesons. The fraction obtained from a fit to the large Mbc 
region is also given. 



Channel 


Fraction, % 


BB 


5.5±J;° 9 ±0.4 


BB*+B*B 


13.7 ±1.3 ±1.1 


B*B* 


37.5 ±\l ±3.0 


Large M bc 


17.5 lij ± 1.3 



To reconstruct the three-body channels, we look for 
an additional charged pion produced directly in the 
3(*)Jj(*) 7r + channels. For each charged pion not in- 
cluded in the reconstructed B candidate, we form right- 
sign B + ir~ , B° ir~ , B° ir + or B~ tt + combinations. We 
then compute the variables M™ s and AE mls for the miss- 
ing B meson, using the energy and momentum of the 
reconstructed Bit combination in the CM: E(B mls ) = 
2E™ m - E(Btt) cm and p(B mis ) = P {Btt) gm . In the 
BB*n+ + B*Bir+ and B* B* ir+ channels, the AE mis 
value will be shifted due to unreconstructed photons from 
B* decays. 

Figure 3(a) shows the corrected Ai? mls + M"" s — m B 
projections for MC simulated BB ir + , BB* n + + B*B n + , 
B*B*ir + , and BBirir events where the B + — > J/ipK + 
mode is generated. The reconstructed B candidates 
are selected from the signal region within the intervals 
5.37GeV/c 2 < M hc < 5.44GeV/c 2 and \AE + M hc - 
m B \ < 0.03 GeV. The value AE mis + M™ is - m B is cor- 
rected by adding to it the value Ai? + Mb c — m B . This 
does not introduce any bias for the original B mesons, but 
improves the resolution, because these two values have 
partially anticorrclatcd uncertainties. Figure 3a shows 
that the BBir+, BB* tt+ + B*Bir+ and B*B* tt+ chan- 
nel contributions are well separated. The reconstruction 
efficiency for the four-body channel (the small peak in 
the rightmost part of Fig. 3(a)) is small and model depen- 
dent. Therefore, we do not include the four-body channel 
in the fit procedure described below. The background 
due to random charged tracks from the unobserved B 
meson is also shown. 

Finally, the A£ mis + A/™ is - m B distribution is ob- 
tained in data for the sum of the five reconstructed B 
modes (Fig. 3(b)). We fit this distribution with a func- 
tion including four terms: three Gaussians with fixed 



shapes and free normalizations to describe the BB ir + , 
BB*tt+ + B*Btt + , and B*B*it+ contributions, and 
a second-order polynomial to describe the background. 
The central positions and widths of the Gaussians are ob- 
tained from fits to the MC simulated distributions shown 
in Fig. 3(a) and are fixed in the fit to data. 
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FIG. 3: (a) The A£ mis + M{£ B - m B distribution normal- 
ized per reconstructed B meson for the MC simulated B + — > 
J/tpK + decays in the (peaks from left to right) BB tv + , 
BB* 7T+ + B*B 7T+, B*B* 7T+, and BB tttv channels, (b) The 
A_B mis ± Afbc is - m B data distribution for right-sign B~ /0 ir + 
combinations for all five studied B modes. The curve shows 
the result of the fit described in the text. 

The event yields and rates are listed in Table III. The 
three-body rates are calculated assuming the ratio of 
charged and neutral directly produced pions to be 2:1 
as expected from isospin conservation. The pion recon- 
struction efficiencies are obtained from the three-body 
phase-space matrix element MC simulation to be 73.8%, 
64.2% and 57.2% for the BBn+, BB* 7t++B*Btt+, and 
B*B* tt + channels, respectively. In neutral B modes, an 
efficiency correction is applied to take into account the 
effect of 19% B a — B° mixing. The systematic uncertain- 
ties due to the fit procedure (±1.5 signal events for each 
channel) and MC efficiency calculations (4-10 %) are also 
included in the total systematic uncertainties. The sta- 
tistical significance of the BB* n + + B*Bir + signal is 
4.4ct. 

It is interesting to note that by directly reconstructing 
pions, we observe only about one half of the rate obtained 
above for the large Mb c region (Table II) resulting in a 
deficit of (9.2 ± 1.0) %. If the difference were due to 
BB 7T7T events, we would expect 4 ± 2 events in the three 
rightmost bins of the AE mls + M™ c ls — m B distribution 
of Fig. 3(b), where only one is observed. In addition, 
the four-body channel is theoretically expected to have 
a rate at least an order of magnitude smaller than the 
three-body channel rates, because there is limited phase- 
space for the creation of an additional pion. Instead, we 
find that the deficit can be explained by ISR contribu- 
tions such as e + e~ — > T(4S)7 — > BB-f. We calculate 
a probability of ~10% for hard photon emission Q by 
the electron or positron beam with subsequent B pro- 
duction, and estimate that ~40 % of such events are due 
to radiative return to the T(4S) resonance. This esti- 
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mate agrees with the observed residual and explains the 
peaking structure on the right side of Fig. 2(b), which is 
dominated by radiative return to the T(4S) or slightly 
higher energies. 

TABLE III: The three-body channel yields and fractions. The 
yields are obtained from a fit to the AE mis + M™ c is - m B dis- 
tribution using five studied B decay modes. The sum of the 
large M^ c channel rates is taken from Table II. The residual 
is the difference between the sum of the three-body channels 
and the result from the large Mb c region. 



Channel 


Yield (tt + ), 


Fraction over 


Fraction 


per 




events 


large M bc % 


bb event 


% 



BBtt 


0.2 


+ 7.2 

-6.9 


r, 2 +6.8 
U ' Z -6.5 


0.0 ±1.2 ±0.3 


BB*tv ± B*Btv 


38.3 


+ 10.5 
- 9.8 


41.6 i£J 


7.3 tH ±0.8 


B*B* 7T 


4.8 


+6.4 
-5.9 


5.9 ±™ 


l.Q±\-£±0A 


Residual 






p 9 q+15.9 
dz -° -15.0 


9.2±l° ± 1.0 


Large M bc 






100. 


17.5 ±i;| ± 1.3 



We analyze other potential sources and backgrounds 
for the events in the multibody region. The rate for a bb 
event to produce the T(4S) and two pions is expected to 
be less than 1 % 0). The wide J p = 1+ B** meson can 
be produced at the T(5S) CM energy with a subsequent 
decay resulting in three- or four-body channels, however 
this process is expected to be negligible due to the very 
small phase-space. The decay 2?° — > B + e~i7 e could con- 
tribute as a background to the multibody channels, how- 
ever, the corresponding branching fraction is estimated 
to be less than 10~ 4 . 

In conclusion, the production of B + and B° mesons is 
measured at the energy of the T(5S). Using fully recon- 
structed B + and B° mesons the production rates per bb 
event are measured to be f(B+) = (72.lj;|;|±5.0)% and 
f(B°) = (77.0±|;| ± 6.1)%. The average value (73.7 ± 
3.2 ± 5.1)% agrees within uncertainties with the CLEO 
value of (58.9 ± 10.0 ± 9.2)% @. Taking into account 
the Bi* ] B^ event rate at the T(5S) of f s = (19.5±"^)% 
fl4| (this value was obtained neglecting bottomonium, 
resulting in an additional absolute uncertainty of about 
1%), and assuming the fraction of final states with a 
bottomonium meson to be 1 - f(B+)/2 - f(B°)/2 - f s , 
some room for unobserved transitions still remains. 

Assuming equal production of B + and B mesons we 
also measure the fractions for bb event transitions to the 
two-body channels with £> + /° meson pairs, f(BB) — 
(5.5 tj;§ ± 0.4) %, f(BB* + B*B) = (13.7 ± 1.3 ± 1.1) %, 
f(B*B*) = (37.5 tfl ± 3.0) %. The BB channel is mea- 
sured for the first time. These fractions are in rough 
agreement with theoretical predictions 0, [l(|, however 
further adjustment of the theoretical models is required. 

Using the additional charged pion directly produced in 
S(*)i?W 7T+ channels, we measure the three-body chan- 
nel fractions in a model-dependent way. The BB* -k ± 



B*Bir decay channel is observed for the first time with 
the fraction f(BB*ir + B*Bir) = (7.3 ±|;f ± 0.8) %. This 
measured three-body fraction is significantly larger than 
those predicted in [ill . [ijj ■ 
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